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Abstract 

Jeans instability is derived for the case of a low density self-gravitating 
gas beyond the Navier-Stokes equations. The Jeans instability criterium 
is shown to depend on a Burnett coefficient if the formalism is taken up to 
fourth order in the wave number. It is also shown that previously known 
viscosity corrections to the Jeans wave-number are enhanced if the full 
fourth order formalism is applied to the stability analysis. 



1 Introduction 

The physics of low density systems is relevant in the description of thermody- 
namic processes which presumably occurred in the universe some time between 
its age of 10 6 sec, temperature T = 1 KeV and the stage after matter decou- 
pled from radiation when the matter density reached values of ~ 10 -18 — and 
T = 1 eV. It is also known that these types of regimes prevail in many other as- 
trophysical systems [Q. Since a hydrodynamical, Navier-Stokes, description of 
a low density fluid breaks down for small Knudsen numbers, it is desirable to go 
beyond this regime and analyze the fluid with the tools that statistical physics 
provide as the particle collision frequency decreases. One simple tool useful for 
this purpose is, on the one hand, the phenomenological formalism now known as 
Linear Irreversible Thermodynamics (LIT), which provides a good grasp about 
how non-equilibrium processes evolve in time |2]. On the other hand, kinetic 
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theory provides a solid mesoscopic foundation for this theory, specially if we 
want to deal with hydrodynamics beyond the Navier-Stokes regime [3] 

The relevance of transport coefficients in cosmological and astrophysical 
problems has already been widely recognized pQ |S] UJ. Nevertheless, the Bur- 
nett equations [3] @] have rarely been applied to structure formation problems 
in a phenomenological context, although the Burnett regime itself has been ad- 
dressed in the context of transport theory in special relativity |H] EH- One 
possible motivation for the study of the Burnett regime in cosmological sit- 
uations is the well-known relation between wave scattering processes, density 
fluctuations and dissipative effects ^U] Dynamic structure factors, 

possibly relevant for describing CMB distortions, can be derived from the hy- 
drodynamical Burnett approach, as it has already been done for collisionlcss 
plasmas The use of these structure factors has appeared to be useful in 

the study of scattering laws to describe the Sunyaev-Zel'dovich effect and can, 
in principle, be applied to study density-density correlation functions to deal 
with hydrodynamical instabilities |14j or temperature-temperature correlation 
functions, a problem that remains to be studied via a thermodynamical ap- 
proach. For the background microwave radiation this is done at present using 
the well-known multipole expansion. In this work we analyze the solutions to 
the phenomenological linearized Burnett equations in the context of Jeans in- 
stability theory, generalizing previous work ^2] JS] . The main reason for doing 
this is that the introduction of the Burnett regime allows to take into account 
all fc 4 order terms in the dispersion relation. Furthermore, the general tenets 
behind this formalism arc worth emphasizing in view of its wide range of appli- 
cability To accomplish this task, the paper is divided as follows: section 
2 reviews the basic phenomenological equations governing the Burnett regime. 
Section three is dedicated to the analysis of the dispersion relation obtained form 
the linearized formalism. Final remarks about the implications of the Burnett 
regime in structure formation are included in the last section. 



2 Basic formalism 

The starting point for our calculation concerns the method whereby for a simple 
fluid which is not in equilibrium we can apply the adequate information to trans- 
form the well-known conservation equations for mass, momentum and energy 
into a complete set. For this purpose, we assume that the fluid is isothermal 
and, therefore, the energy equation may be ignored. This leaves us with the 
continuity and momentum equations. The second one has the form, 

du l dT ij M 

where /* is the external force per unit of volume. To relate the stress tensor r u 
with the mass density p and the local velocity u l , one must appeal to experiment 
or to kinetic theory. The main problem here is that Eq. together with the 
continuity equation involve 4 variables, p and u % . Thus, there are four equations 
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but ten unknowns assuming that the fluid is isotropic and, therefore, r iJ is 
symmetric. So, we require knowledge of the constitutive equations relating t ij 
to the set of local variables. When this procedure is accomplished, keeping all 
those terms which are at most quadratic in the gradients of p and u l we refer 
to it as the Burnett regime. How to accomplish this using kinetic theory is now 
briefly described. We assume that the fluid is a dilute inert gas whose dynamics 
is governed by the Boltzmann equation [3] 0|- For such a model, it is well known 
that through the Chapman-Enskog method one can obtain the solutions to the 
equation which are consistent with the macroscopic hydrodynamic equations. 
This is achieved by expanding the single particle distribution function in power 
series of the Knudsen parameter e, which is a measure of the spatial gradients 
present in the gas, responsible for its deviations from the equilibrium states. 
When this procedure is carried over, one obtains, to order zero in e the Euler 
regime, and to first order in e the Navier-Stokes regime with zero bulk viscosity, 
so that 

4 s =p6v -2^ --66*) (2) 

where <r y is the symmetric traceless part of the velocity gradient, r\ is the shear 
viscosity, p is the local pressure, 5^ is Kronecker's delta and 9 = To 
next order in e one gets the Burnett regime, in which t ij becomes a rather 
complicated expression pi] containing both, nonlinear as well as linear terms 
in the velocity gradient. If we further assume that the gas is slightly deviating 
form its local equilibrium state, the linearized Burnett contribution to the stress 
tensor Tg is then given by: 

2 d(-^-f) i d dv 1 d{—P) 1 
Po oxj p oxj axi 3 axj p a 

In this equation n72 ia a transport coefficient whose value is known for rigid 
spheres and Maxwell molecules as well, p and p are the equilibrium values for 
the density and pressure, respectively and p the local pressure p(x l ,t). t v in 
Eq. is the total stress tensor, T z ^ s + Tg , where t^ s is given in Eq. J2J. When 
we add these contributions, take their divergence and use the local equilibrium 
assumption, p = p(p, T = const), so that = |^ J^- — C 2 S ^~ one finally gets 
that the linearized Burnett approximation to Eq.Q is given by: 

- iv.v ■ (I) - A c ;v 2 v,(^)i = /' (4) 

Po Po 3 po Op 

where C 2 S is the speed of sound. 

The linearization procedure, which is identical to the one followed in a pre- 
vious publication |14j . has clearly invoked that local state variables p and u 1 
deviate form its equilibrium values through linear fluctuations, namely 

p = p + Sp (5) 
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and 

u l = u* + 5u l = 6u l (6) 
The continuity equation now reads as 

d(Sp) , 



at 



(7) 



where 69 = Eqs. Q-Q constitute the set of linearized coupled equations 
for bp and u % which can be solved for a given conservative force /*. In our case 
of interest, the external force density is taken to be the gravitational force and, 
moreover, we assume that the fluctuating part of the gravitational field Sip is 
governed by Poisson's equation 

V 2 5cp = AnGSp (8) 
where G = 6.67 x lO -11 /^ is the gravitational constant. 

kgs z o 

To eliminate u % from this set we simply take the divergence of Eq. 
and make use of Eq. © plus the fact that |V' = V • -Jj so that after a 
straightforward arrangement of terms we find that 



3 p \ dt J 3 p po 

Eq. © is the main result in this calculation. It is a single equation describing 
the density fluctuations of an isothermal fluid when subject to the action of 
a fluctuating gravitational field taking into account all linearized contributions 
up to second order in the Knudsen parameter. Its implications on the Jeans 
number will be studied in the following section. Before we do so, however, it 
is interesting to point out the main features that in Eq. J^J originate from the 
Burnett approximation, those proportional to rj 2 . The one of fourth order in the 
spatial derivatives, the next to the last term in the left hand side, corresponds 
to what was to be expected. But the other one namely, the third term in the left 
hand side shows a coupling between the gravitational field and the dynamics of 
the density fluctuations. This coupling is totally absent both at the Euler and 
the Navier-Stokes level. Its implications will be discussed later. 



3 Analysis of the dispersion relation 

In order to solve Eq. @ we propose that the density fluctuations are represented 
by plane waves, so that 

5p(r,t) =Ae- i{wt -^ (10) 

where the amplitude A needs not to be specified and the rest of the symbols 
have their standard meaning. Substitution of Eq. Ill 01) into Eq. yields: 
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uj 2 + k 2 (-C 2 + ^W2V 2 G Po ) + l^iuok 2 - ^ 2 ^/fc 4 + 4nGp = (11) 
3p 3 p 3 popo 

Eg. llljl is the sought dispersion relation from which we must examine the con- 
ditions leading to values of ui, necessarily complex, giving rise to runaway terms 
in Eq. I)l()|l . In fact, it may be rewritten as 



where 



and 



uj 2 + 2ia(k)u + f(k) = (12) 

a{k) = \^k 2 (13) 
3 p 

C 2 8tt 2 2 C 2 ru 2 r] 2 A 

f(k) = (— ~ q — Gw2w) - - — -k 4 + ATiG Po (14) 

7 3p Q 3 7p D p 

When examining the roots of this quadratic equation we see that the runaway 
solutions require that uj be complex and the threshold values for k to meet this 
requirement is that a(k) 2 + f(k) > 0. This condition reduces to a quadratic 
equation for the minimum value of k 2 namely, 

r 2 k 4 + sk 2 +4nGp = (15) 

where 

r = q — - Gp (16) 

y Po 6 

S = -C 2 -^ 2V 2 ^- (17) 
3p 7Po 

Calculating the roots of Eq. I|15fl . taking the root with the positive sign, 
k 2 > 0, and expanding the square root we find, after a few trivial algebraic 
steps, that 

TT = H k J\ 7777^3 — ) ( 18 ) 

kj po JK Po 9Cs 2 Po J 

where the ordinary Jeans wave number is given by kj = ^2o. Eq.JHJ is the 
sought result. In its derivation we have consistently neglected all terms of order 
rf and higher. 

To what extent this result has an influence on the ordinary Jeans number will 
be discussed below. Two minor points should be noticed. Firstly, if voi = the 
corrective term, which is of order of fc 4 is the one obtained within the Navier- 
Stokes approximation when ignoring the bulk viscosity |16j . Secondly, as we 
have pointed out initially, the Burnett linearized terms have to be included to 
have a consistent correction to order fc 4 . 



and 
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4 Discussion 



As Eq. i|18|) clearly indicates, the inclusion of the Burnett linearized term pro- 
vides a contribution to the correction of fcj which is opposite in sign to the one 
arising form the Navier-Stokes equations. This may have a significant effect on 
the Jeans mass, which is inversely proportional to kj [Sj ^2] |15| . 

One may now ask about the conditions in which dissipative effects are indeed 
relevant in structure formation problems. In order to address this question we 
start with typical monatomic gases such as H, He, Ar, etc which we know, are 
reasonably described by a hard-sphere model whose shear viscosity extracted 
from kinetic theory is given by 0]: 



ksmT 



" = is?V^r- (19) 

Here a is the atomic diameter, m the mass of the particles, ks is Boltzmann's 
constant and T is the temperature of the system. The most relevant contribution 
will arise from the term within brackets. To estimate its value we recall that 
for hard spheres zu 2 — 2 and the equation of state is p = P °^ T , where R is 
the universal gas constant and M the molecular weight of the gas. Then, the 
correction is determined by the expression 

" - ^T ( ~^r ~ vcjio* (20) 

Using the equation of state, the value of w% and the fact that for a monatomic 
gas Cf = | , then for hydrogen M = 2 whence Eq. PU)| reduces to 

Since 7 ~ 1.7 for monatomic gases, this equation finally reads as 

8 a 4 R Po T Po T v ' 

Eq. I122II has to be taken as a result which requires some care before con- 
clusions may be extracted from it. As we said in the introduction, the Burnett 
regime is an adequate formalism to deal with fluids in the so called transition 
regime. This implies the regime prevailing between the rarefied regime and the 
continuous one. In the latter, the well-known Navier-Stokes equations of hy- 
drodynamics provide an accurate description of the fluid. It appears that this 
transition regime could have prevailed in the Universe for temperatures T ~ 10 4 
(K) and densities 10~ 15 -Ss, but these figures are merely qualitative. If such 
conditions existed, then Eq. I|22|l simply asserts that structure formation would 
have been enhanced due to the coupling between hydrodynamic and gravita- 
tional modes. Notice that this coupling is the dominant term of Eq. (|21() . 
Clearly the importance of the correction critically depends on p D and T. If p Q T 
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is larger than 10 -19 as it probably occurred in earlier stages, the correction C 
will be negligible so that, in order to reach a more definite statement, more 
reliable data is required. As a final comment, it should be mentioned here that 
in this calculation we have ignored taking into account the expansion in a co- 
moving reference frame. Whether or not this improves the results is a subject 
of further study. 

This work has been partially supported by CONACyT (Mexico), project 
41081-F. 
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